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ABSTRACT 

By means of chemo-photometric models for galaxies of different morhological, we have 
carried out a detailed study of the history of element production by spheroidal and 
dwarf irregular galaxies. Spheroidal galaxies suffer a strong and intense star formation 
episode at early times. In dwarf irregulars, the SFR proceeds at a low regime but 
continuously. Both galactic types enrich the IGM with metals, by means of galactic 
winds. We have assumed that the galaxy number density is fixed and normalized 
to the value of the optical luminosity function observed in the local universe. Our 
models allow us to investigate in detail how the metal fractions locked up in spheroid 
and dwarf irregular stars, in the ISM and ejected into the IGM have changed with 
cosmic time. By relaxing the instantaneous recycling approximation and taking into 
account stellar lifetimes, for the first time we have studied the evolution of the chemical 
abundance ratios in the IGM and compared our predictions with a set of observations 
by various authors. Our results indicate that the bulk of the IGM enrichment is due 
to spheroids, with dwarf irregular galaxies playing a negligible role. Our predictions 
grossly account for the [O/H] observed in the IGM at high redshift, but overestimate 
the [C/H]. Furthermore, it appears hard to reproduce the abundance ratios observed 
in the high-redshift IGM. Some possible explanations are discussed in the text. This 
is the first attempt to study the abundance ratios in the IGM by means of detailed 
chemical evolution models which take into account the stellar lifetimes. Numerical 
simulations adopting our chemical evolution prescriptions could be useful to improve 
our understanding of the IGM chemical enrichment. 

Key words: Galaxies: evolution; Galaxies: intergalactic medium; Galaxies: funda- 
mental parameters. 



1 INTRODUCTION 

In the universe, metals can be found in different amounts in 
various environments. A fraction of all the metals produced 
so far is located in galaxies, locked up into living stars and 
stellar remnants, or contained in the inter-stellar gas. The 
remaining fraction of metals has been ejected by galaxies 
into the intergalactic and intracluster media (IGM, ICM). 
The amount of metals present in galaxies and IGM/ICM 
is a function of the cosmic epoch. From the high-redshift 
universe to the present time, the observations indicate an 
evolution of the metallicity of galaxies (e.g Kobulnicky & 
Kewley 2004). In lockstep with the growth of the galactic 
metallicity, the universe has experienced a progressive metal 
enrichment of the inter galactic gas. Metal abundances in 

* E-mail: fcalura@ts.astro.it 



the IGM are measurable by means of the quasar absorption 
lines (see Bechtold 2003, Pettini 2004). At high redshift, the 
observations report IGM metal abundances between lO"'^'^ 
and 0.1 solar (Lu et al. 1998, Cowie & Songaila 1998, Berg- 
eron et al. 2002, Carswell 2004, Simcoe et al. 2002, Schaye 
et al. 2003, Simcoe et al. 2004a, Aguirre et al. 2004). In the 
local universe, the available data indicate an average IGM 
metallicity of the order of ~ 0.1 solar (Tripp et al. 2002, 
Shull et al. 2003). 

Many theoretical studies have attempted to model the global 
chemical evolution of the universe and the chemical evolu- 
tion of the IGM (Pei & Fall 1995, Edmunds & Phillips 1997, 
Malaney & Chaboyer 1996, Pei, Fall & Hauser 1999, Gen 
& Ostriker 1999, Sadat Guiderdoni & Silk 2001, Mathlin 
et al. 2001, De Lucia et al. 2003, Dunne et al. 2003), very 
fruitful in tracing the average properties of the universe but 
unsuited to study in detail the different contributions to the 
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overall metal content of the universe from the single galactic 
types, namely spheroids, spiral discs and irregulars. Further- 
more, many galactic evolution models compute metal pro- 
duction by adopting the instantaneous recycling approxima- 
tion (IRA), which docs not allow one to study the elements 
whose production occurs on timescales not negligible rela- 
tive to the age of the universe, such as C and Fe. In a recent 
paper, by means of chemical evolution models for gala^cies 
of different morphological types, Calura & Matteucci (2004, 
hereinafter CM04) have computed the total amount of met- 
als produced by ellipticals, spirals and irregular galaxies and 
locked up in all the different components of the present-day 
xmivcrsc. 

In the present paper, by means of the same detailed galac- 
tic chemical evolution models, we investigate how the metal 
fractions in the IGM have changed with cosmic time. One 
of the main eissumptions of CM04 was that all galaxy types 
started forming at high redshift but with quite different star 
formation histories. In particular, elliptical galaxies arc as- 
sumed to have formed in a very short timescale (less than 1 
Gyr) and therefore, since they are the major metal produc- 
ers, they are the most likely candidates to enrich the IGM 
at high redshift. A large catlicgory of objects, experiencing 
very intense starbursts are observed at high redshift, thus 
confirming our assumption. These objects, which include the 
SCUBA (Blain et al. 1999, Trentham, Blain & Goldader 
1999) and ultra luminous infrared galaxies (ULIRG, Rowan- 
Robinson 2000), very luminous in the infrared band but faint 
in the optical, dominate the cosmic star formation at red- 
shift z s (Perez-Gonzalez et al. 2005) and are natural 
candidates for massive proto-ellipticals (Ivison et al. 2000). 
On the other hand, our assumptions are at variance with the 
hierarchical clustering scenario for galaxy formation. One of 
the main predictions of this scenario is that small struc- 
tures, intended as dark matter halos, form first. As a con- 
sequence, the first galaxies to form are the smallest ones. 
The largest galaxies arc formed later on by means of several 
merging episodes among gas rich galaxies, such as discs, oc- 
curring through the whole Hubble time (Kauffmann, White 
& Guiderdoni 1993, Baugh et al. 1998, Cole et al. 2000, 
Somerville et al. 2001, Menci et al. 2002). These models 
were supported by some observational results like the ap- 
parent paucity of massive galaxies at 2 ~ 1 — 2 (Barger et 
al. 1999, Menanteau et al. 1999, Fontana et al. 1999, Drory 
et al. 2001, Dickinson et al. 2003). However, recent results 
seem to show that there is a consistent number of ellipti- 
cals at those redshifts (Schado et al. 1999, Im et al. 2002, 
Yamada et al. 2005). Beside spheroidal galaxies, we study 
also the contribution that dwarf irregular galaxies brought 
to the IGM metal enrichment. 

We try to give a possible answer to several fundamental 
questions in galactic chemical evolution, such as: when was 
the IGM enriched with metals? Can massive spheroids form- 
ing at high redshift account for the metals in the IGM? Is 
a "normal" (i.e. Salpeter-like) initial mass function (IMF) 
sufficient to enrich the IGM at the observed values? 
This paper is organized as follows: in section 2 we recall all 
the main assumptions of the chemical evolution models and 
the method used to compute the comoving metal densities. 
In Section 3 we present the results and in Section 4 we draw 
the conclusions. 



2 THE CHEMICAL EVOLUTION MODELS 

We assume that the main contributors of the IGM enrich- 
ment are elliptical galaxies, the main spheroidal component 
of spiral galaxies, i.e. bulges, followed by dwarf irregular 
galaxies. Elliptical galaxies, spiral bulges and halos are all 
included in the cathegory of spheroids, and in this work are 
described by one single model. 

In this paper, we assume that spiral discs do not develop 
galactic winds and do not contribute to the chemical enrich- 
ment of the IGM. This assumption is motivated by several 
reasons. First, in the disc of our own Galaxy, winds do not 
seem to have played an important role in its evolution. In 
fact, models of chemical evolution of the galactic disc with 
mass outflows do reproduce the abundances and abundance 
ratios observed in field stars (Tosi et al. 1998, Matteucci 
2001). In general, the kincmatical features of the interstellar 
gas in spiral discs seem rather different from the ones ob- 
served in local and distant objects experiencing winds and 
outflows. Zaritsky et al. (1990) have measured the veloci- 
ties of Hii in local spirals, finding values of ay ~ lOfcms^^. 
These values are much lower than the ISM speeds in local 
starburst and Lyman-Break Galaxies, in which very ener- 
getic material moving at speeds > 250fems~^ has been de- 
tected (Pettini et al. 2002). The gas velocities measured in 
local spirals do not seem thus able to generate violent phe- 
nomena such as strong outfiows or winds. 
We calculate galaxy magnitudes by means of the spectro- 
photometric code by Jimenez et al. (1998) This code allows 
us to compute self-consistently the photometric evolution of 
composite stellar populations taking into account the metal- 
licity of the gas out of which the stars form at any time (see 
CM03 for a more detailed description). 



2.1 The spheroidal galaxy model 

The chemical evolution model representing spheroids (ellip- 
ticals and bulges) has been developed in order to reproduce 

a large set of observablos concerning elliptical galaxies, such 
as abundance ratios and photometric properties (Matteucci 
1994, Pipino & Matteucci 2004) . Detailed descriptions of the 
chemical evolution models used in this work can be found 
in Matteucci (1992, 1994). 

In our picture, spheroids form as a result of the rapid col- 
lapse of a homogeneous sphere of primordial gas where star 
formation is taking place at the same time as the collapse 
proceeds. Star formation is assumed to halt as the thermal 
energy Eth of the ISM, heated by stellar winds and super- 
nova (SN) explosions, balances the binding energy of the 
gas. At this time a galactic wind occurs, sweeping away al- 
most all of the residual gas. By means of the galactic wind, 
ellipticals enrich the IGM with metals. 
The thermal energy Eth of the ISM is calculated according 
to: 

Eth = {Eth)sN + {Eth)w (1) 

[EthJsN represents the energy injected by SNe and is given 
by: 

{Eth)SN = esNIa^^O^isNIa + eSNII^o J?SNII (2) 
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where Eq = 10^^ erg is the energy released by each SN, esNia 
(esNii) and 7?sNia (Rsmi) are the efficiency of energy transfer 
and the explosion rate for SNe la (SNe II) , respectively. Our 
model takes into account also radiative cooling and follows 
the prescriptions suggested by Cioffi, McKoo & Bcrtschinger 
(1988) for the SN remnants evolution (see Pipino et al. 2002, 
2005). According to this formulation, the net amount of en- 
ergy transferred from SNe into the ISM is not larger than 
20 % of the initial blast wave energy. {Eth)w indicates the 
energy injected by stellar winds and is given by: 



(j){m)tp{t )r]wEwdmdt 



(3) 



Stellar winds give a negligible contribution to the total ther- 
mal energy of the gas. Ew ~ 3-10*'^ erg is the energy released 

into the ISM from a typical massive star (~ 2OM0) through 
stellar winds during its lifetime and rjw is the efficiency of 
energy transfer, which we assume to be 3% according to 
Bradamante et al. (1998). 

(p(m) represents the stellar initial mass function (IMF), 
which has been assumed to be constant in space and time. 
Unless otherwise stated, we assume a Salpeter IMF, which 
is given by: 



0(m) = m' 



i-l+x) 



(4) 



where 0o is a normalization constant and x = 1.35 over 
all the stellar mass range, mi„/ < m/MQ < ruup, with 



O.IMq and niu 



IOOMq. We will test also the 



possibility of an IMF flatter than the Salpeter one in the 
low-mass end, similar to the one found in the local disc of 
the MW (Chabrier 2003, Kroupa 2002). 
ip{t) is the star formation rate (SFR), namely the fractional 
amount of gas turning into stars per unit time. The SFR 
ip{t) (in Gyr~^) is given by: 



ip{t) = vGit) 



(5) 



The quantity G{t) = Mg(t)/Mtot is the total fractional 
mass of gas present in the galaxy at time t, being Mg(t) 
and Mtot{t) the gas mass and total meiss at the time t, re- 
spectively. The quantity v is the efficiency of star formation, 
namely the inverse of the time scale for star formation and 
for a typical elliptical with 10^^ M0 of baryonic matter is 
assumed to be ~ 10 Gyr~^. As shown by Matteucci (1994), 
to reproduce the observed correlation between the magne- 
sium index Mg2 and the iron index < Fe > in ellipticals, 
this quantity has to be higher in more massive galaxies. This 
leads to the so-called "inverse wind" scenario, in the sense 
that large galaxies develop a wind before the small ones. 
The SFR is assumed to drop to zero at the onset of a galac- 
tic wind. 

The binding energy Esit) is given by the expression (see 
Pipino et al. 2002): 



EBit) = WLit) + WDMit) 



(6) 



where Wl (t) represents the gravitational energy of the gas 
due to the luminous matter: 



WL{t) = - 



0.5GMg{t)ML{t) 



(7) 



where Mg{t), ML{t) and rz, represent, at the instant t, the 
gas mass, the total luminous mass and the effective radius. 



respectively. WDM{t) is the gravitational energy due to the 
dark matter (DM), and is given by: 



WdmH) 



GMg(t)MDark{t) 



Wr 



(8) 



Muarkit) = WMl is the amount of DM present in the 
galaxy, W^m takes into account the distribution of the dark 
matter relative to the luminous one: 

w'nM{t) = ^^^[\ + \.2.7{rL/ru)] (9) 

Vh/rD is the ratio between the effective radius and the core 
radius of the DM. We set Th/rD = 0.1, in agreement with 
the results by Bertin et al. (1992). 

The prescriptions for the DM halo adopted here, namely 
a rather heavy but diffuse DM halo, are supported by sev- 
eral observational results (Gerhard et al. 2001, Romanowsky 
et al. 2003, Samurovic & Danziger 2005), which have shown 
that in several early type galaxies, there is no strong ev- 
idence for dark haloes out to several effective radii. The 
adopted dark matter profile is hence different than a cuspy 
profile, like the one suggested by Navarro, Frenk & White 
(NFW, 1996). Several recent theoretical results seem to ex- 
clude a cuspy DM profile or presenting a peak in the center 
as the one by NFW (Mamon & Lokas 2005, Mcrritt et al. 
2005). It is worth to note that in SPH simulations of el- 
liptical galaxies, galactic winds in very massive galaxies are 
difficult to explain owing to deep dynamical potential in the 
central regions (Kobayashi 2005). It is possible that in such 
simulations, the adoption of a more diffuse DM profile, like 
the one of this paper, could solve this problem. Concerning 
the baryonic matter, we take into account the self gravity of 
the gas and we assume that the gas distribution is the same 
as the one for the stars at any time. Let Gi be the fractional 
mass of the element i in the gas within a spheroidal galaxy, 
its temporal evolution is described by the basic equation: 



Gi = -iji{t)Xi{t) + Riit) - {Gi)out 



(10) 



If Mg{t) is the mass of the gas at the time t, Gi{t) = 
Mg{t)Xi{t)/Mtot is the gas mass in the form of an ele- 
ment i normalized to the total initial mass Mtot- If G(t) = 
Mg{t)/Mtot, the quantity X^{t) = Gi{t)/G{t) represents the 
abundance in mass of an element i, with the sum over all 
the elements in the gas mixture being equal to unity. Ri{i) 
represents the returned fraction of matter (both newly pro- 
cessed and already present in the star) in the form of an 
clement i that the stars eject into the ISM through stellar 
winds and SN explosions. This term contains all the pre- 
scriptions regarding the stellar yields and the SN progenitor 
models, and is given by: 



Ri{t) = Sml"' ~ '^rn)Qmi{t - rm)(t>{m)dm+ 

Mb. 

Mu 



+ (1 - A) J^^^ i;{t - Tm)Qmi{t - Tm)<p{m)dm+ 



I 



i>{t - Tm)Qmi{t - Tm)(j){m)dm 



(11) 



Each one of these integrals represents the contribution 
of different stellar mass ranges, including type la SNe (sec- 
ond integral). In particular, Qmi(t — Tm) = Qmi X{t — Tm), 
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where Qmi is a matrix which calculates for any star of a 
given mass m the amount of the newly processed and al- 
ready present element i, which is returned to the ISM. The 
quantity Tm is the lifetime of a star of mass m. The constant 
A is the fraction of binary systems which produce type la 
SNe. For all the details of the calculations of the type la 
SN rate, we address the reader to Greggio & Renzini (1983) 
and Matteucci & Greggio (1986). The term {Gi)out accounts 
for the presence of gas outflows, occurring by means of SN- 
drivcn galactic winds. {Gi)out is simply equal to the rate 
of mciss loss from dying stars after the occurrence of the 
wind, since all the matter present in the ISM is lost on a 
timcscalc of a few 10* years. We assume that metals ejected 
into the IGM by any galactic type do not accrete on the 
other galactic types. Recent studies have indicated that, if 
the star formation is rapid (< 0.5 Gyr), the presence of an 
infall term has negligible effects on the predicted chemical 
abundances (Pipino & Matteucci 2004). For this reason, we 
assume that all the gas is already present at the beginning 
of star formation, which is equivalent to a very short infall 
episode. 



2.2 The dwarf irregular galcixy model 

The dwarf irregular galaxies are assumed to assemble all of 
their gas by means of a continuous infall of pristine gas. 
Hence, the chemical evolution equation is slightly different 
than eq. 10, since it contains also an infall term, which is 
given by: 



{Gi)inf — C 



{Xi)iufe 



<t/T) 



Ml 



(12) 



where {Xi)i„f is the abundance of the element i in the in- 
falling gas, assumed to be primordial, r is the time scale of 
mass accretion and C is a constant obtained by imposing to 
reproduce the total luminous mass Ml at the present time 
to. The parameter r has been assumed to be the same for 
all dwarf irregulars and short enough to avoid unlikely high 
infall rates at the present time (r = 0.5 • l(f years). 
As in the case of ellipticals, dwarf irregulars are assumed 
to be embedded in a dark matter halo. The prescriptions 
for the dark matter halo are similar to those adopted for 
the ellipticals, namely a massive but diffuse dark halo (see 
Bradamante et al. 1998 for details). The SFR proceeds at a 
low regime but continuously. The SFR expression is identi- 
cal to the one assumed for ellipticals (see eq. 5 ), but the SF 
efficiency adopted in this case is z/ = 0.05Gyr~^ . The SF, if 
strong enough, can also trigger a galactic wind, in principle 
easy to develop in irregular galaxies owing to their shallow 
gravitational potential wells. The conditions for the onset 
of the wind are identical to the ones described for elliptical 
galaxies (Bradamante et al. 1998, Recchi et al. 2002). 
Also in this case, we assume a Salpeter IMF. The chem- 
ical evolution model for dwarf irregulars allows us to re- 
produce the main features of local blue compact dwarfs 
(BCD) and dwarf irregular galaxies. The predicted present- 
day ISM O abundance for the typical dwarf galaxy model 
is 12+log(0/H)= 8.5, with local observations of BCDs in- 
dicating values between 7.15 and 9 (Shi et al. 2005). The 
predicted HI mass to blue luminosity ratio is Mhi/Lb = 
O.5M0 /Lq , with the observations indicating values of 0.2 — 



2Mq/Lq (Garland et al. 2004). The predicted SFR at the 
present time is O.lMo/yr, while local observation in irreg- 
ular magellanic and BCDs show values of 0.01 — IMo/j/r 
(Hopkins et al. 2002, Hunter & Elmegreen 2004). 



2.3 Stellar Yields 

In both the models used in this work, the chemical enrich- 
ment of the ISM occurs by means of three processes: stellar 
winds, type la and type II SN explosions. Low and inter- 
mediate (0.8 < M/Mq < 8) mass stars enrich the ISM by 

means of quiescent mass loss, or stellar winds. For these kind 
of stars, we adopt the nucleosynthetic prescriptions by van 
den Hoeck & Groenewegen (1997). 

Type la SN are assumed to occur in binary systems with 
total masses 3 < Mb/Mq < 16, formed by a white dwarf 
(WD) accreting mass from a secondary star until the Clian- 
drasekhax mass (~ 1.4Mo) is reached (Greggio & Renzini 
1983). At this point, the star explodes by means of the ther- 
monuclear runaway and no remnant is left. 
Single massive (M > 8AIq) stars enrich the ISM by means 
of type II SN explosions. After the explosion, a neutron star 
or a black hole can be left as remnants. For massive stars 
and type la SNe, we adopt the empirical yields suggested by 
Frangois et al. (2004), which are substantially based on the 
Woosley & Weaver (1995) and Iwamoto et al. (1999) yields, 
respectively. For C production, we do take into account 
also rotation in massive stars, consistently with Meynet and 
Maeder (2005). Therefore, the results for these elements are 
slightly different than the ones obtained in Calura & Mat- 
teucci 2004. 



2.4 Comoving densities of chemical elements and 
IGM abundances 

To study galaxy evolution, we assume a pure luminosity evo- 
lution picture. According to this assumption, galaxies evolve 
only in luminosity and their number density is constant in 
time. Such a picture can account for many observables, such 
as the evolution of the galaxy luminosity density in various 
bands and the cosmic supernova rates (Calura & Matteucci 
2003, hereinafter CMOS). 

The galactic morphological types considered in this work 
are two: spheroids and dwarf irregulars. The production rate 
Pi{z), for the i—th element and the galaxies of A;— th mor- 
phological type as a function of the redshift z is given by: 



pUMq yr-' Mp<r^){z) = p),{z){M j L)l{zH {z) 



(13) 



where p\ {£) is the B band luminosity density generated by 
galaxies of fc— th morphological type and is the result of the 
integral over all luminosities of the morphological B-band 
luminosity function (LF) ^^{Lb)(z), namely: 



{Lb){z) {LB/L*B{z))dLB 



(14) 



The general LF #(Ls) is parametrized according to the form 

defined by Schechter (1976): 

^Lb){z) = {Lb/L*b{z))-'' exp{-LB/L'B{z)) (15) 
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where $* is a normalization constant related to the number 
of luminous galaxies per unit volume, Lg (z) is a character- 
istic luminosity and a is the faint-end slope. 
At redshift zero, we adopt the morphological LF as observed 
by Marzkc ct al. (1998). At redshift larger than zero, wc have 
calculated the LF by applying the evolutionary corrections 
calculated by means of the spectrophotometric model by 
Jimenez et al. (1998) to the local LF. 

The evolution of the stellar, gaseous and total mass is a di- 
rect prediction of the chemical evolution models, whereas 
the evolution of the B luminosity Ls is calculated by means 
of the spectro-photometric code. (M/L)%{z) represents the 
B mass-to-light ratio for fc— th galaxies, where M stands for 
the total baryonic mass, (z) represents the rate of pro- 
duction of the i— th element by a stellar generation (only 
the newly processed fraction) in fe— th galax;ies. At a given 
time t, it is given by: 



rl' 

^ m 



ipk{t - Tm) Tmi{t - Tm) ^(rn)dm 



(16) 



The quantity Tmi is the matrix Qmi previously defined, with 
the diagonal elements (material already present in the star) 
set to zero. The conversion between the cosmic time t and 
redshift z is performed on the basis of the adopted cosmolog- 
ical model, which in our case is characterized by Q,m = 0.3, 
Q.A = 0.7 and Ho = 65Km s"^ Mpc'^. 
At a given redshift z, the total comoving mass density in the 
form of a given element (labeled i) produced in the galaxies 
of A;— th morphological type ^^((MoMpc^'^) is given by 
the integral: 



Pi,Tot{z) 



• k / f 
Pi (z 



dt 
dz' 



dz 



(17) 



where p'i{z') is the production rate for the i-th element as a 
function of redshift for the galaxies of fc— th morphological 
type. 

As in CM04, the metals produced by spheroids and 
dwarf irregulars end up into three phases: some are locked 
up in stars plus remnants, some arc contained in the ISM 
gas and the rest is expelled into the IGM. At each instant, 
we can predict the comoving density of metals contained in 
each phase. To calculate the comoving mass density for the 
element i and for fe— th galactic type in stars, ISM gas and 
IGM, first we need to calculate the average stellar and ISM 
abundances predicted by our models for a typical spheroidal 
and dwarf irregular galaxy. These quantities are direct out- 
puts of our chemical evolution models. 

Let us define Xi ,,{z) and Xi^g{z) as the average abun- 
dances by mass locked in stars and remnants and in ISM 
gas, respectively, as a function of the redshift z for the ele- 
ment i in A;— th galactic type. The stellar comoving density 
of the element i in fc— th galactic type is given by: 



pX*{z) = Xl,{z) ■ pl{z) 



(18) 



whereas the comoving density of the element i in the ISM 
gas is given by: 



piiSM{z) = Xlg{z) ■ pI{z) 



p^{z) and p^{z) are the mass densities of stars plus remnants 
and ISM in the galaxies of the fc— th morphological type, 
respectively, and are calculated as: 

p':{z)=pI{z)-{M./Lb)\z) (20) 
and 

p1{z)=p'h{z)-(MJLB)\z), (21) 

where {M^/Lb)''{z) and {Mg/LB)''{z) represent the pre- 
dicted stellar and gas mass to light ratios for fc— th galactic 
type, respectively. 

The comoving density of a given element i expelled by 
the galaxies of the fc— th morphological type into the IGM 
is given by: 



Jo 



pIigm{z) = pUz)-{MJLbT{z)- I iGi)tt{z')^dz' (22) 



The total mass density in the IGM for the element i is 
then given by: 



pi,iaM{z) = y^^PijGMjz) 

k 



(23) 



We assume that, immediately after the enrichment, the met- 
als are distributed in the IGM homogeneously, i.e. we assume 
that the IGM is "well mixed" . This is probably a rough ap- 
proximation, since the physical scale is very large, of IMpc 
or more. Furthermore, not only metal enrichment, but also 
gas cooling and metal diffusion depend on the environment, 
i.e., local densities. These effects are not taken into account 
in our model, but we compensate by calculating the chemi- 
cal enrichment in great detail. 

To determine pb,iGM we have proceeded as in CM04. The 
first thing we need is a value for the total baryonic density. 



= phi pc 



(24) 



where pc = 1.17 • 10^^ Mq/Mpc^ is the critical density of 
the universe, calculated for the value of /i = 0.65 adopted 
here. For the baryonic density, here we assume the fiducial 
value of Qt — 0.02h~^ (Spergcl et al. 2003), corresponding 
to pB ^ 5.55 • 10'-^ Me/Mpc^. To calculate the IGM baryon 
density, we then subtract to this value the baryon density in 
galaxies (stars and ISM) as calculated in CM04. The bary- 
onic density in galaxies can be constrained by means of the 
observed luminosity density, as shown by Calura, Matteucci 
& Menci (2004). They have shown that the assumption that 
the bulk of ellipticals formed at high redshift allows us to 
reproduce a large set of observational data. 
We determine the abundances in the IGM in the following 
way: 



lXi/H]iGM{z) = log 



Pi.IGMjz) 

Ph,igm{z) 



- log{Xi/H), 



(25) 



(19) 



where pi,iGM{z) is the mass density of the element i which 
has been ejected in the IGM by spheroids and dwarf irregu- 
lars and Ph,igm{z) is the comoving density of H in the IGM. 
To determine ph,igm{z), for the IGM we assume a primor- 
dial composition with ^i^e abundance of 0.25 (see Steigman 
2005) and H abundance of 0.75. 

{Xi/H)Q is the solar abundance of the element i relative to 
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Ht. 

In the present work we do not consider metal production by 
pop III stars, and we assume that only normal stars in high- 
redshift galaxies are the sources of the metals residing in 
the IGM. This assumption is motivated by the results from 
Matteucci & Calura (2005), who have shown that pop III 
stars enrichment cannot account for the abundance ratios 
observed in the high redshift IGM. Furthermore, the critical 
metallicity necessary to enter a regime of normal star forma- 
tion is Zcr = (10"*' - 10"'')^Q (Ferrara & Salvaterra 2004, 
Yoshida, Bromm, Hernquist 2004), i.e. well below the values 
observed in the IGM considered here, namely Z > W~^Zq. 

3 RESULTS 

3.1 The evolution of the cosmic metal share 

In Figures 1 and 2 we show the predicted metal comov- 

ing densities contained in the stars, in the ISM and those 
ejected into the IGM as functions of redshift, as produced 
by spheroids and dwarf irregulars, respectively. In Figures 
3 and 4 we show the evolution of the metal fractions in the 
ISM, IGM and in the stars, i.e., at each redshift, the ratio 
between the metal mass density in each phase (stars, IGM 
and ISM) and the total one for any given element, as pro- 
duced by spheroids and dwarf irregulars, respectively. Here, 
we show only the results obtained by assuming that both 
spheroids and dwarf irregulars started forming stars at red- 
shift Zf = 5. Figure 1 and 2 provide information on the 
values of the total comoving mass densities, expressed in 
MqMpc~^ , obtained by summing up the various phases at 
any redshift. On the other hand. Figures 3 and 4 are useful 
to see how the fractions of metals contained in the various 
phases change with the redshift. 

In the case of galactic spheroids (Figures 1 and 3), we see 
that elements produced on different timescales show differ- 
ent behaviours. Elements produced by low and intermediate 
mass stars, such as C and N, are mainly produced on long 
timescales (> IGyr), with their bulk being produced after 
the end of star formation. 

Elements produced by massive stars (M > 8Mq) on short 
timescales (< O.OSGyr), such as O and Mg, are produced 
during the star formation epoch. At the present time, nearly 
half of these elements is locked up into the stars. 
Fe is an elements produced by both type la and II SN explo- 
sions, but in different proportions: 3/4 of the Fe originates 
in type la SNe, whereas the remaining 1/4 comes from type 
II SNe (Matteucci & Grcggio 1986). According to the nu- 
cleosynthesis prescriptions adopted here, Zn is produced in 
a substantial way in SNe la. Therefore, both Fe and Zn are 
then massively produced after the starburst. 
Si is an element produced both by type II and type la SNe 
and in similar proportions. At the present time, nearly 30 % 
of all the Si produced in spheroids is locked up into stars. 
In the case of dwarf irregular galaxies, the total metal den- 
sity and stellar metal density (solid and dotted lines in Fig- 
ure 2, respectively) is a slowly increasing function through- 
out all their history. This is due to a continuous and low 

t The solar abundances used in this paper are taken by Anders 
& Grevesse (1989). 



star formation rate characterizing these objects. The ISM 
metal densities (dot-dashed lines in Figure 2) increase until 
the time at which the galactic wind develops, i.e. for the 
model used here, 3.6 Gyr after the beginning of star forma- 
tion, corresponding to redshift z — 1.5. At this redshift, the 
metal densities in the ISM start to decrease, whereas the 
metal densities in the IGM start to increase. 

3.2 Chemical abundances in the IGM 

The abundances in the IGM are determined by means of the 
metal absorption lines observed in the quasar (QSO) spec- 
tra. In these spectra, the presence of the IGM is indicated 
by the multitude of individual Lyman-a (Lya) absorption 
lines bluewards of the Lya emission peak of the QSO. These 
absorption lines originate in ionized gas clouds with typical 
column densities in the range 10^** > N{HI) > 10^^ cm"^ 
(Pettini 2004, Bechtold 2003), located in the hues of sight 
of the QSOs. The abundances in the Lya forest are usu- 
ally determined for C, O and Si. Many are the uncertain- 
ties in the determination of the abundances and concern 
the modeling of the IGM conditions (Simcoe et al. 2004b, 
Carswell 2004). Major uncertainties are also the nature of 
the UV photon flux, i.e. whether visible or dust-obscured 
QSOs (Rees & Setti 1970, Pei & Fall 1995) or star forming 
galaxies (Haardt & Madau 1996, Kim, Cristiani, D'Odorico 
2002, Bianchi, Cristiani & Kim 2001) represent the principal 
sources of UV photons. Observationally, the of the Oxygen 
abundance [0/H] is derived by using measurements of Ovi 
and Hi absorptions, according to: 

[O/H],,, = log{^) + log{^) - log{^)Q (26) 

where iVo^j and Nhi are the column densities of the chem- 
ical species Ovi and Hi, respectively, /o^/ and /hj are the 
ionization fractions, i.e. the ratios between the column den- 
sity of the ionized species Ovi and Hi with respect to the 
neutral ones (Simcoe et al. 2004a). The same method may 
be applied to calculate [C/H] and [Si/H] from the Civ and 
Siiv absorption lines, respectively. The main challenge in 
determining the abundances is an accurate estimate of the 
ionization fractions, depending on all the physical properties 
of the IGM discussed above. 

In figures 5 and 6 we show the evolution of the intergalactic 
C (upper panel) and O (lower panel) abundances as pre- 
dicted by our models and as observed by various authors, 
having assumed that all galaxies started to form stars at 
Zf = 3 and Zf — 5, respectively. For this comparison we 
consider only these two elements as they represent the ones 
with the largest number of measures. Some of the observa- 
tions considered here have been derived by testing differ- 
ent shapes of the ionizing UV backgrounds (dominated by 
QSO or a mixture between galaxies and QSOs). These data 
arc represented in Figures 5 and 6 by squared and shaded 
regions. Data derived by assuming one single type of UV 
background are represented by points with error bars. 

In these figures we show both the contributions by 
spheroids (thin solid lines) and dwarf irregulars (thin dotted 
lines) to the total IGM abundances, along with the sum of 
the two (thick dashed lines). For both elements, it is worth 
to note that the contribution by dwarf irregular galaxies to 
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Figure 1. Predicted redshift (lower x-axis) and lookback-time (upper x-axis) evolution of the metals produced by spheroids (various 
elements). Solid lines: total mass density pi produced in tlic form of a given clement i (here the letter Z indicates the total mctallicity). 
Dotted lines; metal mass density locked up in stars. Dot-dashed lines: metal mass density contained in the ISM. Dashed lines: metal 
mass density ejected into the IGM. In this case, we assume that in all spheroids star formation starts at redshift Zf = 5. 
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Figure 2. Predicted redshift (lower x-axis) and lookback-time (upper x-axis) evolution of the metals produced by dwarf irregulars 
(various elements). Solid lines: total mass density pi produced in the form of a given clement i (here the letter Z indicates the total 
metallicity ) . Dotted lines: metal mass density locked up in stars. Dot-dashed lines: metal mass density contained in the ISM. Dashed 
lines: metal mass density ejected into the IGM. In this case, we assume that in all dwarf irregulars star formation starts at redshift 
Zf = 5. 
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Figure 3. Predicted redshift (lower x-axis) and loolcbaclc-time (upper x-axis) evolution of the metal fractions, normahzed to the total 

mass in the form of a given clement, locked up in galactic spheroids in stars (dotted lines), in the ISM (dot-dashed lines), in the IGM 
(dashed lines) for various elements. In this case, we assume that in all spheroids star formation starts at redshift Zf = 5. 



the IGM appears negligible with respect to spheroids. Dwarf 

irregulars contribute to the 6% of the C mass density and 
to the 12% of the O mass density present in the IGM at the 
prsent time. This is an important result, which is in agree- 
ment with the ones already found in Calura & Matteucci 



(2003), where we showed that these objects have played a 
negligible role in the cosmic star formation history. In the 
case of [C/H], in both figures 5 and 6 we note that at high 
redshift the predicted total IGM abundances are in agree- 
ment with the value by Bergeron et al. (2002, light-hatched 
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Figure 4. Predicted redshift (lower x-axis) and lookback-time (upper x-axis) evolution of the metal fractions, normalized to the total 
mass in tlic form of a given clement, locked up in dwarf irregular galaxies in stars (dotted lines), in the ISM (dot-dashed lines), in the 
IGM (dashed lines) for various elements. In this case, we assume that in all dwarf irregulars star formation starts at redshift Zj = 5. 



area), but overestimate the IGM [C/H] abundance as mea- 
sured by Simcoc ct al. (2004a, square box) and Schaye ct 
al. (2003, dotted area). If we take the observational data at 
face value and neglect all their uncertainties, we note that 
our predictions calculated by assuming Zf = Z seem in better 



agreement with the measures than the ones calculated by as- 
suming Zj — 5. The predictions are in very good agreement 
with the estimates by Barlow & Tytler (1998, shaded area) 
at redshift z ~ 0.5, but underestimate the value by Tripp 
et al. (2002, empty diamond) at redshift zero by ~ 0.4 dex. 
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Figure 5. Evolution of the C (upper panel) and O (lower panel) abundances in the IGM as measured by various authors and as predicted 
by our models by assuming that the star formation in spheroids and dwarf irregulars starts at Zf = 3. Thin solid lines: contribution by 
spheroids. Thin dotted lines: contributions by dwarf irregulars. Thick dashed lines: total abundances. In both panels, the observational 
data indicated with squared and shaded regions have been derived by taking into account different ionizing UV backgroud shapes. Data 
in upper panel: open diamond: Tripp et al. (2002). Solid triangle; Barlow & Tytler (1998). Light-hatched area: Bergeron et al. (2002). 
Square box: Simcoe et al. (2004a). Dotted area: Schaye et al. (2003). Data in lower panel: open diamond: Tripp et al. (2002). Open 
hexagon: ShuU et al. (2003). Solid square: Simcoe et al. 2002. Empty square: Dave et al. (1998). Dashed area: Telfer et al. (2002). Crossed 
box: Carswell et al. (2002). Square box: Simcoe et al. (2004a). 



The one by Tripp et al. (2002) is the only local observational 
estimate of the C abundance in the IGM, and it has been 
derived by considering only a UV background dominated 
by QSOs (Haardt & Madau 96). A softer UV background, 
including also galaxies, would produce significantly lower C 
abundances (see Schaye et al. 2003). 

The most reliable IGM metallicity tracer is represented by O 
since, at variance with C, it is observable down to the lowest 
density regions of the universe, characterized by HI column 
density values of Nhi < W^"^'^ cm~^ . Unfortunately, also the 
determination of the abundances of O suffers the same un- 
certainties as C. Furthermore, O is detected by means of the 
Ovi absorptions, which arise in the Lyman a forest. For this 
reason, the measure of its column density is contaminated 
also by blending effects. 

The [O/H] in the IGM (figures 5 and 6, lower panels) 
is roughly accounted for by our predictions both at high 
and low redshift. We note that the observations show no 
strong evolution in the O abundances between z = 2.5 and 
z = 0, suggesting that the IGM was substantially enriched 
with metals at very high redshift. At high redshift {z > 3), 
our predictions are consistent with the value by Dave et al. 
(1998, empty square), which has been provided with no er- 
ror bar. However, if the uncertainties are the same as for 
the data by the other authors, i.e. S > 0.5 dex, we can say 
that this measure is more consistent with our predictions 



calculated by assuming a redshift of formation zf — 5. Our 
predictions overestimate the values observed by Carswell et 
al. (2002 crossed regions in figures 5 and 6) and Simcoe et 
al. (2004a, large square box), but are in very good agree- 
ment with the data by Telfer et al. (2002, dashed region). 
The only measured abundance being slightly higher than our 
predictions is the value by Simcoe et al. (2002) at z ~ 2.5. 
Although being underestimated by our models and higher 
than any other value, this measure represents a lower limit 
and it is unlikely to represent the average IGM O abun- 
dance. It is likely that this object represents a galactic wind 
in action, or a proto-cluster region: in fact, the intergalac- 
tic systems associated to those absorptions are located in 
highly overdense environments, where the gas density is par- 
ticularly large and the star formation and metal production 
rates could be enhanced. 

Finally, the data at redshift zero by ShuU et al. (2003, open 
hexagon) and the lower limit by Tripp et al. (2002) are in 
very good agreement with our predictions. 
If the main responsible of the IGM enrichment are the big 
spheroids, this lack of evolution in the IGM O abundance 
is another independent indication that the largest spheroids 
have formed at very high redshift, consistently with many 
other observational evidences such as the fundamental plane 
(Renzini & Ciotti 1993), the Mg-(j» relation (see Bernardi 
et al. 2003 and references therein) and the color-magnitude 
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Figure 6. Evolution of the C (upper panel) and O (lower panel) abundances in the IGM as measured by various authors and as predicted 
by our models by assuming that the star formation in spheroids and dwarf irregulars starts at z:^ = 5. Lines and data as in Figure 5. 



relation (Bower, Lucey & Ellis 1992). 

From the analysis described in this section, we conclude that 
our predictions can account for the observed evolution of the 
[0/H] in the IGM, being consistent with most of the exis- 
tent data. In the case of [C/H] in the IGM, our predictions 
overestimate 2 of the 3 measures quoted above. This could 
be due to inadequate ionization corrections applied to the 
Civ column density estimates. It is well known that the C 
abundances are based primarily on Civ observations, which 
probe the densest 10% of the baryons in the universe and the 
1% of the cosmic volume (Simcoe et al. 2004b). In regions 
where the density is near the universal mean. Civ loses its 
sensitivity as a metallicity tracer, since most of the C has 
been ionized to Cv and higher levels. In those regions, we 
could be still missing the bulk of the C. If this is true, all 
the observational data should be regarded as lower limits. 
On the other hand, although more difficult to detect owing 
to blending effects, in these very low density regions Ovi 
does not lose its reliability as a metallicity tracer, owing to 
its ionization potential and the high O abundance (Carswell 
et al. 2002). 

Another possible reasons for the discrepancies between pre- 
dictions and data for both C and O could be the fact that 
the IGM is not well mixed as assumed in the present work. 
The data could be biased and could not represent the aver- 
age values of the whole universe. Finally, a high fraction of 
the metals could not b observable, residing at higher tem- 
perature in the intra-cluster hot gas, or in warm gas. 



3.3 The IMF and the abundance ratios in the 
IGM 

The possibility of a non-standard Salpeter IMF in various 
environments has already been suggested by various authors. 
In the disc of the Milky Way, the stellar IMF is steeper than 
the one typical of young open clusters (Kroupa 2004; Scalo 
2004). This is supported by results from chemical evolution 
models, finding that in general a Salpeter-like IMF tends 
to overestimate the metal abundances observed in the local 
disc (Romano et al. 2004; Portinari et al. 2004). In elliptical 
galaxies, it has been noted that the standard Salpeter IMF 
tends to overestimate the average M/Lb and that the best 
results are obtained by flattening the IMF at masses < IMq 
(Renzini 2004). 

From a theoretical point of view, the main parameters affect- 
ing the abundance ratios are the adopted stellar yields and 
the stellar IMF. Here, we assume a fixed set of yields and 
we try to modify the IMF, in order to investigate also the 
effects of a non-standard Salpeter-like IMF on the predicted 
abundance ratios. In figure 7, we compare observational es- 
timates of the [Si/C] (upper panel) and [0/C] (lower panel) 
with our predictions, by assuming that star formation in 
spheroids started at redshift Zf = 5. We are showing predic- 
tions calculated assuming various types of IMF: a standard 
Salpeter (solid lines) and a modified Salpeter with flatter 
slope, X = 0.35 instead of a; = 1.35, below two mass-cuts, 
O.6M0 (dotted lines) and 2.OM0 (dashed lines), respectively. 
Recently, Aguirre et al. (2004) have carried out a detailed 
study of silicon, in the IGM in the redshift range « ~ 2 — 4. 
By using the results presented by Schaye et al. (2003) , who 
measure carbon abundance in the same redshift range, and 
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assuming different UV background types, they have pro- 
vided an observational estimate of the [Si/C] ratio in the 
IGM. The [Si/C] observed at high redshift by Aguirre et al. 
(2004, shaded area in Figure 7, upper panel) spans a broad 
range, from [Si/C]=0.25, obtained with a UV background 
(UVB) including gala:x;ies, to [Si/C]=1.5 in the case of UVB 
only due to QSOs. As a redshift range, we have assumed 
the one of the QSOs studied by Aguirre ct al. (2004), i.e. 
1.5 < 2 < 4.5. A local estimate of the [Si/C] ratio in the 
IGM is provided by Tripp ct al. (2002), who studied the 
chemical abundances in two Lya clouds in the vicinity of 
the Virgo supercluster, finding [Si/C] ~ 0.2 ± 1 (empty di- 
amond in upper panel of figure 7) . 

Our predictions arc consistent with the range of values pro- 
vided by Aguirre et al. (2004). We note that no choice of 
the IMF allows us to have values for the [Si/C] ratio at high 
redshift larger than [Si/C] ~ 0.4 dex. 

At redshift zero, the predicted [Si/C] values for various IMFs 
arc higher than the value observed by Tripp et al. (2002). 
More observations with a detailed study of the effect of dif- 
ferent ionizing UV background shapes would be very useful 
in order to understand the nature of this discrepancy. 
In Figure 7, lower panel, we show the evolution of the [O / C] 
ratio in the IGM, as predicted by our models and as observed 
by various authors. The data for the [O/C] ratio in the IGM 
at high redshift {z > 1.5) are from Bergeron et al. (2002, 
Light-hatched area) and Telfer et al. (2002, shaded area). At 
redshift zero, we use the values by Tripp et al. (2002, empty 
diamond), where [0/H] represents a lower limit. At high 
redshift, our predictions lie outside of the range indicated 
by Telfer et al. (2002), but are marginally consistent with 
the data by Bergeron et al. (2002). In general, the flatter is 
the IMF in the low mass regime, the higher is the obtained 
[O/C] value and the better is the agreement between obser- 
vations and predictions. Globally, our results indicate that it 
is hard to reproduce the abundance ratios [O/C] and [Si/C] 
observed in the IGM at rcdshifts larger than 2. One way 
to match the observations would be to adopt a top-heavy 
IMF in spheroids, but in CM04 we have already shown that 
such a choice would lead to a severe overestimation of the 
local metal budget. On the other hand, the adopted stellar 
yields could be inappropriate. The predicted carbon abun- 
dance is too large. In order to reproduce the observed [Si/C] 
and [O/C], it would be necessary to decrease ad- hoc the C 
yields adopted in this work. However, it is worth to stress 
that for C, the nucleosynthesis yields adopted here allow us 
to reproduce the abundance pattern observed in the Galac- 
tic stars with good accuracy (Chiappini et al. 2003). With 
such an ad-hoc modification of the C yields, it could be hard 
to reproduce the local C abundance. 

Another possible solution could be that a non-negligible 
fraction of C is locked up into dust grains. This would 
cause an underestimation of the actual IGM C abundance. 
The treatment of dust formation in the chemical evolution 
equations could be very useful to investigate this possibility 
(Calura et al., in preparation). Finally, inhomogencitics and 
local effects could play an important role, not taken into 
account in our models. 



3.4 Escape of metals in IGM and stall radius 

An important aspect of the IGM metal enrichment concerns 
how far the metals can travel in the space once they are 
ejected, and whether they can reach the lowest density re- 
gions of the IGM where metals are detected. Ferraxa, Pettini 
& Shchekinov (2000) have performed a theoretical study of 
the mixing of the metals in the early IGM. By means of a 
semi-analytical model for SN feedback, they concluded that 
SN-driven winds, or blowouts, fail to distribute the met- 
als over volumes large enough to pollute the whole IGM. 
Both dwarf and giant galaxies eject metals over maximum 
distances Rmax of < IQOkpc, a factor of 10 smaller than the 
average physical separation between massive objects at high 
redshifts, which is of the order of IMpc. This means that we 
should expect that regions outside the radius Rmax have a 
primordial chemical composition. Ferrara et al. (2000) con- 
clude that SN-driven galactic winds arc inadequate to ex- 
plain the ubiquitous presence of metals in the IGM, and 
that additional mechanisms must be at work to account for 
the efficient mixing of the metals on scales of ~ IMpc. 
On the other hand, in a more recent paper, Tumlinson & 
Fang (2005) have studied the distribution of the intergalac- 
tic Ovi absorbers in the local IGM. They concluded that, 
in order to account for the observed distribution of metals 
in the local IGM, massive galaxies have to enrich regions of 
~ 0.5 - IMpc. 

On the basis of simple physical arguments, we attempt to 
calculate the maximum distance, or "stall radius " Rs which 
could by reached by the ejected metals. 
We assume that a spheroid of baryonic mass Mlbg ~ 
lO^^M© and Mulirg = 10^^ Mq represents a typical Lyman 
break galaxy (LBG) and ULIRG, respectively. According to 
the inverse wind scenario, we assume that the star formation 
efficiency increases with the galactic mass: for the 10^° M© 
and 10^^ A/0 galaxies, we assume star formation efficiencies 
I'LBG = 3 Gyr^^ and vulirg = 11 Gyr^^ (Matteucci 1994, 
Pipino & Matteucci 2004). 

In Figure 8, we show the SFR, expressed in M^/yr^ for 
the models representing the ULIRG (solid line) and the 
LBG (dotted line), as a function of time. According to our 
results, a 10^^ M© elliptical develops a wind in 0.23 Gyr, 
with average values for the SFR of ~ 1O"M0/2.3 • 10*j/r ~ 
400M©/j/r. The SFR of the most massive elliptical galaxies 
can reach very high peak values (up to ~ 1000MQj/r~^), 
similar to the ones observed in some SCUBA (Ivison et al. 
2000) and luminous infrared galaxies (Rowan-Robinson et 
al. 1997, 2000). 

On the other hand, a lO^^Mo spheroid develops a wind 
in ~ 0.6 Gyr, hence the average star formation rate is 
~ VJMQ/yr. 

Here, the basic assumptions for the ISM heating by SNe 
and stellar winds are the same eis described in various pa- 
pers (Pipino et al. 2002, Pipino et al. 2005). As in Pipino et 
al. (2002), we assume that every SN releases an energy = 
10 '^ erg, ~ 20% of which goes into thermal heating of the 
ISM, after taking into account SN remnant cooling. We as- 
sume that all the thermal energy is carried out by the wind. 
In this case, at each timestep the energy lost in the galactic 
wind is Ai? = Eth, where Et^ is given by equation (1). If 
the quantity AM represents the amount of matter ejected 
through the galactic wind at each timestep, we can then cal- 
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Figure 7. Evolution of the [Si/C] and [O/C] abundance ratios in the IGM as predicted by our models for different IMFs and as measured 
by various authors. Upper panel: solid line: Salpeter IMF. Dotted line: IMF flatter (x = 0.35) than the Salpeter at masses m < O.GMq. 
Dashed line: IMF flatter {x = 0.35) than the Salpeter at masses m < 2Mq. Observations: shaded area: Aguirre et al. (2004); open 
diamond: Tripp et al. (2002). Lower Panel: curves as in upper panel. Observations: shaded area: Teller et al. (2002). Light-hatched area: 
Bergeron et al. (2002). Open diamond: Tripp et al. (2002). In this case, we assume that star formation in spheroids starts at 2:^^ = 5. 
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culate the wind velocity according to — {2AE / AM)^^^ 
(Pipino et al. 2005). With the prescriptions adopted here, 
immediately after the starburst we obtain wind velocities of 
Vw,LBG ^ 500Km/s and Vw,ulirg ~ lAOOKm/s in the case 
of a typical LBG and ULIRG, respectively. Both values are 
compatible with the ones measured in real LBGs (Pettini et 
al. 2002, Lehmer et al. 2005), local starbursts (Martin 2005) 
and LIRGs/ULIRGs (Veilleux 2003, Veilleux et al. 2005). 
To give an estimate of the stall radius Rs , we assume that the 
wind propagates as a spherical bubble in a uniform medium. 
We require that the outflow ram pressure balances the IGM 
pressure pi cm (Ferrara et al. 2000), obtaining: 

Rs=(^^Y^ (27) 

where is the mass loss rate of the galaxy. According to 
our predictions, immediately after the onset of the wind the 
mass loss rate is ~ 1.9M0/t/r and QM^/yr in the case of 
the LBG and ULIRG, respectively. The pressure of the IGM 
is given by: 

p/GM =p* [(! + «) /200]'' (28) 

where pt/Kn = 1500 h^cm~^K is the IGM pressure at 
z = 200 (Ferrara et al. 2000) . A fundamental quantity to cal- 
culate piGM is the IGM temperature Tigm- The structure of 
the IGM is likely to be multiphase, with a colder medium at 
temperatures T = lO"^-^'^ K (Danforth & ShuU 2005, Stocke, 
ShuU & Penton 2004). Here, we assume for Tigm the value 
in-between these two media, i.e. Tigm = lO^K. With these 
hypothesis, at rcdshift z = 4 for the stall radius we obtain 
a value for Rs,lbg = 0.77Mpc and Rs,ulirg = 2.8Mpc 
for LBGs and ULIRGs, respectively. We conclude that, on 
the basis of these very simple arguments, the metals ejected 
by massive galaxies can reach very long distances, being 
in principle able to account for the observed ubiquity of 
the metals in the diffuse IGM. Furthermore, it is worth 
to note that, with such values for the wind velocity and 
assuming this quantity to be constant, the ejected mat- 
ter would be able to move to the required distances of 
~ IMpc on timescales of ~ lMpc/vw,LBG = 1.95Gyr and 
^ IMpc/v w, ULIRG — 0.7 Gyr J in the case they are ejected 
by a LBG or a ULIRG, respectively. The most massive galax- 
ies, experiencing the most intense starburst in the universe, 
could be responsible for the enrichment of the lowest density 
regions of IGM, which are also the most distant from galax- 
ies. It has been found recently that the role of ULIRGs in the 
cosmic star formation increases rapidly at redshift z > 1.5, 
and that at rcdshift z ~ 3 they completely dominate the 
SFR density (Perez-Gonzalez et al. 2005). Furthermore, the 
presence of very massive starbursts is detected up to z ^ 5.5 
(Staguhn et al. 2005), along with the presence of very metal- 
rich massive quasars (Dietrich et al. 2003, Iwamuro et al. 
2004). On the light of these evidences, it is likely that these 
objects dominate also the metal enrichment of the IGM at 
redshift z ~ 3 — 5. In principle, by assuming that very in- 
tense starbursts have occurred at high redshift , it would 
then be possible to explain the presence of metals in very 
undcrdcnsc regions of the IGM even at rcdshift z~ 3 — 5, 
i.e. when the universe was only ~ 1.5 Gyr old. We want to 
stress that these conclusions have been reached on the ba- 
sis of pure order-of-magnitude arguments, and that only de- 



tailed numerical simulations following in detail the transport 
of metals in the IGM could give more precise estimates of 
the quantities described above. It would be very interesting 
to probe whether, with the prescriptions adopted above, in 
particular star formation rates of the order of > 400AIq /yr, 
speeding up the ejecta at very high velocities, such simula- 
tions could be able to account for the presence of the metals 
on scales of ~ IMpc or more. 



4 CONCLUSIONS 

By means of chcmo-photomctric models for spheroidal and 
dwarf irregular galaxies, we have carried out a detailed study 
of the history of element production in the universe. Accord- 
ing to these models, the galaxy densities are fixed and are 
normalized to the values of the optical luminosity function 
in the local universe (Marzke et al. 1998). This is equiva- 
lent to assume that galaxies evolve in luminosity and not in 
number. In our picture, the metal enrichment of the IGM is 
mainly due to galactic spheroids (i.e. ellipticals and spiral 
bulges) and, to a lesser extent, to dwaxf irregular galaxies. 
We did not consider galactic winds in spiral discs. Our main 
results can be summarized as follows. 

i) The enrichment of the IGM has proceeded in different 

manners for elements produced on different timescales. In 
spheroids, a-elements (O, Mg) arc produced during the star- 
burst phase, whereas for elements synthesized by low and 
intermediate mass stars (C, N, Fe) the bulk of the produc- 
tion occurs well after the starburst. The behaviour of Fe, 
produced mainly by type la SNc, is similar to the one of the 
elements produced by low and intermediate mass stars. Zn is 
also perhaps produced by type la SNc, but a non-negligible 
contribution comes by massive stars (Frangois et al. 2004). 
Dwarf irregular galaxies, characterized by a prolonged star 
formation history, in general have a different metal share 
than spheroids, with a low fraction locked up in stars (less 
than - 10%), ~ 20% in the ISM and the highest fraction ( 
~ 70%) ejected into the IGM. 

ii) The bulk of the IGM enrichment is due to spheroids, with 
dwarf irregular galaxies playing a negligible role. Our predic- 
tion can reproduce the [O /H] observed in the IGM by various 
authors, as well as the trend of [O/H] with redshift. On the 
other hand, the predicted [C/H] are larger than observed, 
especially at high-redshift, in 2 out of 3 cases. A possible 
explanation could be that in regions of the universe with 
density close to the cosmic mean, where the vast majority 
of the baryons lie. Civ loses its sensitivity as a metallicity 
tracer, since C has been ionized to higher levels, whereas 
Ovi is a rrmch more sensitive metallicity probe. Other pos- 
sible reasons for the discrepancies between predictions and 
data for both G and O could be the fact that the IGM is not 
well mixed, or the data could be biased and representative 
of the average cosmic mean, or a a high fraction of the met- 
als could not be observable, residing at higher temperature 
in the intra-cluster hot gas, or warm gas. 

iii) The study of the abundance ratios may be useful to con- 
strain the IMF of the first stars in galaxies. We have tried to 
fit the [Si/C] and [0/C] ratio observed in the local and high- 
redshift IGM by assuming a standard Salpeter IMF and a 
modified Salpeter flattened below various possible mass val- 
ues. Our predictions underestimate the [0/C] and [Si/C] 
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Figure 8. Star formation rate, expressed in Mq/uv, for an elliptical galaxy of baryonic mass M = IO^^Mq (solid line) and M = IO^^Mq 
(dotted line), as a function of time. 



observed in the IGM at redshift larger than 2. To reduce 
this discrepancy, the adoption of a top-heavy IMF is not a 
viable solution, since it would overestimate the metal bud- 
get in the local universe (CM04). The explanation resides 
perhaps in inhomogeneities and local effects, which could 
be investigated by means of numerical simulations studying 
single elements and taking into account the stellar lifetimes, 
iv) On the basis of simple but realistic physical arguments 
we have shown that, in order to disperse the metals up to 
distance scales of ~ IMpc in less than IGyr, outflow veloci- 
ties of vw ~ 1500Km/s for the galactic winds are required. 
These velocities are of the same order of magnitude of the 
values observed in luminous infrared galaxies, experiencing 
the most intense starbursts in the universe. If present also 
at high redshift, as many recent observations seem to sug- 
gest (e.g. Perez-Gonzalez et al. 2005), these objects could be 
responsible for the enrichment of the lowest density regions 
of the IGM. 

We stress that at the present time, the observations of the 
metal abundance in the IGM have been obtained for a few 
elements (mainly C and O). By improving the number of 
detections for these elements and by deriving abundances 
for other ones, it will be possible to shed more light on the 
sources of the IGM metals and on the mechanisms governing 
the IGM enrichment. 
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